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A numerical solution of Wheeler-De Witt equation for a quantum cosmological model 
simulating boson and fermion creation in the early Universe evolution is presented. This 
solution is based on a Wheeler-De Witt equation obtained by Krechet, Fil’chenkov, and 
Shikin, in the framework of quantum geometrodynamics for a Bianchi-I metric. 
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It is generally aserted that in the early stage of Universe evo- 
lution, the quantum phase predominated the era. Therefore 
there are numerous solutions have been found corresponding 
to the Wheeler-DeWitt equation which governs this phase [2]. 
In the present paper we present another numerical solution of 
Wheeler-De Witt equation for a quantum cosmological model 
simulating boson and fermion creation in the early Universe 
evolution for a Bianchi-type I metric [1]. 

The solution is based on Wheeler-De Witt equation for a 
Bianchi-I metric obtained by Krechet, Fil’chenkov, and 
Shikin [1], in the framework of quantum geometrodynamics. 
Albeit the essence of the solution is quite similar from the so- 
lution given in [1] using Bessel function, in the present paper 
we present numerical result using Maxima. For comparison 
with other solutions of 1-d hydrogen problem, see [3] and [4]. 


2 Solution of Wheeler-DeWitt equation for boson and 
fermion creation 


In the evolution of the Universe after inflation, a scalar field 
describing de Sitter vacuum was supposed to decay and its 
energy is converted into the energy of fermions and heavy 
vector-particles (the so-called X and Y bosons) [2]. 
In the framework of quantum geometrodynamics, and for 
a Bianchi-I metric, the Wheeler-De Witt equation has been 
obtained by Krechet, Fil’chenkov, and Shikin, which reduces 
to become (Eq. 23 in [1]): 
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(E-V)T =0. (1) 
where T” and T’ represent second and first differentiation of 
T with respect to r. The resulting equation appears quite 
similar to radial 1-dimensional Schrödinger equation for a 
hydrogen-like atom [3], with the potential energy is given 


by [1]: B ” 
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has here a continuous spectrum. 

The solution of equation (1) has been presented in [1] 
based on modified Bessel function. Its interpretation is that 
in this quantum cosmological model an initial singularity is 
absent. 

As an alternative to the method presented in [1], the nu- 
merical solution can be found using Maxima software pack- 
age, as follows. All solutions are given in terms of Æ as con- 
stant described by (3). 


(a) Condition where V = 0 
diff(y,r,2) -E*y—(2* %i*C/3/t)*y=0; ode2(%ol,y,r); (4) 


The result is given by: 





y=K1 sin (a)+K2 cos (a), (5) 
where: 
a=(r/ V3) /—3E—2iC/t. (6) 
(b) Condition where V $ 0 
'diff(y,r,2)—E*y—(2*%i*C/3/t)*y—(b/t+e/t*/3)*y=0; 
ode2(%02,y,r); (7) 
The result is given by: 
y=K1 sin (d)+K2 cos (d), (8) 
where: 
d=(r/(/3 2/2) V—3Ett/3—2iCt!/3—3e—3bt!/3. (9) 





As a result, the solution given above looks a bit different 
compared to the solution obtained in [1] based on the modi- 
fied Bessel function. 
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3 A few implications 


For the purpose of stimulating further discussions, a few im- 
plications of the above solution of Wheeler-DeWitt equation 
(in the form of 1-d Schrödinger equation) are pointed as fol- 
lows: 


(a) Considering that the Schrödinger equation can be used 
to solve the Casimir effect (see for instance Silva [5], 
Alvarez & Mazzitelli [6]), therefore one may expect 
that there exists some effects of Casimir effect in cos- 
mological scale, in a sense that perhaps quite similar to 
Unruh radiation which can be derived from the Casimir 
effective temperature. Interestingly, Anosov [7] has 
pointed out a plausible deep link between Casimir ef- 
fect and the fine structure constant by virtue of the en- 
tropy of coin-tossing problem. However apparently he 
did not mention yet another plausible link between the 
Casimir effective temperature and other phenomena at 
cosmological scale; 


(b) Other implication may be related to the Earth scale ef- 
fects, considering the fact that Schrédinger equation 
corresponds to the infinite dimensional Hilbert space. 
In other words one may expect some effects with re- 
spect to Earth eigen oscillation spectrum, which is re- 
lated to the Earth’s inner core interior. This is part of 
gravitational geophysical effects, as discussed by Gr- 
ishchuk et al. [8]. Furthermore, this effect may corre- 
spond to the so-called Love numbers. Other phenom- 
ena related to variation to gravitational field is caused 
by the Earth inner core oscillation, which yields oscil- 
lation period T ~ 3-7 hours. Interestingly, a recent re- 
port by Cahill [9] based on the Optical fibre gravita- 
tional wave detector gave result which suggests oscil- 
lation period of around Shours. Cahill concluded that 
this observed variation can be attributed to Dynami- 
cal 3-space. Nonetheless, the Figure 6c in [9] may be 
attributed to Earth inner core oscillation instead. Of 
course, further experiment can be done to verify which 
interpretation is more consistent. 
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